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I. INTRODUCTION 

The spin-1/2 XXZ chain is one of the fundamental 
models in the study of the low-dimensional magnetism. 
The Hamiltonian is given by 

N 

ff = ^{5J5jVi + ^j'^jVi + A5|5|+i}, (1) 

where S — a/2 and a are Pauli matrices. The model 
can be solved by Bethe ansatz method and diverse 
physical properties have been investigated with vary- 
ing anisotropy parameter A P, |3j I3- For A = 1 we 
know the second neighbor correlation 0|. Jimbo et al 
y, Q obtained the multiple integral representation for 
the arbitrary correlation functions of the XXZ chain for 
1 > A > — 1 and A > 1. Later Kitanine, Maillet and 
Terras ,5] generalized to magnetized case using the Bethe 
wave function. But this multiple integral is very com- 
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plicated and they only reproduce the nearest neighbor 
correlations. Recently Boos et al succeeded to calculate 
these multiple integrals for emptiness formation proba- 
bihty (EFP) 

m ^ 

P(m) = (n(^|+2))' (2) 

at A l,m = 3,4,5,6[i[ll[ll[l3. The case to = 3 
reproduces the result of the second neighbor correlation. 
The same method is applied to the third neighbor cor- 
relations at A = 10. In the previous paper we calcu- 
lated the second neighbor corelation functions for mass- 
less case — 1 < A < 1 and integral is simplified to the one 
dimensional integral [l^. In this paper we give the sec- 
ond neighbor correlation in massive region 1 < A. The 
correlation formula contains the elliptic functions. But 
we can show that the final result also contain only one- 
dimensional integral. 
■i |siroisi@issp. u- tokyo . ac . j p| 
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II. MULTIPLE INTEGRAL FORMULA 



A. Massless case 



The correlation function in massless case is given by, 

j./^i,^2,...,.^v .... . ._^'r-^^dx, - r dx^^. sinh^(A.-A.) 

^e'l' 4, -,4^ (V'glV's) ^ C' /J y-oo-.c 27r ^.^ii ^ 27r J-^i sinh(A, - A, - 

yr smy-\n'j + ?,^^)s\n\r-'{^i'^ + i^) ^ sinh^-^(A^j - 4) sinh'"-^(/^j + z§) 
^ -^J- (coshf^')™ -^J- (coshf/i,)™ ' 



where 



A = cosC, {Al, A„J = {^j,^^^, nj,^^^^, fij^^^, fij^^J, 

a+ = {j; 1 < j < m, e,j = +}, card(Q!+) = s', MaXjg„+ (j) = J^aa;- Minjg„+ (j) = 

a" = {j; 1 < j < TO, e'j = -}, card(a") = s, Max^g^- (j) = jmaa:, Min^g^- (j) = j^m. (4) 



If we put Xj Xj — for j = 1, s', these equations are written as follows 

i^D':?'-':r)=(-ir(-J)^ 

veil ^21 tm^ C 



c 

n sinh 5(A„ - A5) n sinh^^"'(A, + sinh'"-^' (A, - 



-00 27r 7rAj_ J]^ sinh(Aa - A^ - ifabC) 



a>b 



where 



fab = (1 + sign[(s' + 1/2 - a)(s' + 1/2 - 6)])/2, 

2/1 > 2/2, ... > 2/s', Cy. = +, Vs' + l > ys'+2 > ... > ym, Em+l-y, = 

If we put Xj —^ XjC, this is written as 



[| sinh''^"-! C{xj + ^) sinh"-^^- C{xj * ' 



a>6 

The weight function of the integral is common for all the correlation functions: 



n(-7r) sinh7r(a;Q - Xb) 

W^{xi, ...,Xm) = ^ = nSech(7ra;j)) = Wi{xj)), 



J]^ 2 cosh™ TTXj 



where Dm is antisymmetric differential operator defined by 



Especially in the limit C — > or A = 1 we have 



l[ix, + ^)y^~\x,~^r-y^ 



j^l-J-oo [[{Xa - Xb - ifab) 



a>b 

These integrals are done analytically up to m = 4 for A = 1 and up to m = 3 for — 1 < A < Ipl 



B. Massive case 

For massive case integral formula for correlation function contains elliptic functions!^ 0,0 



\e„e^,...,e„J LLJ ^^^^^^^ 2n ^^J-^J _^/^ 2n t>b''^'^^^^ 



Xb 



X 



-pj- sin^-i(/i;- +3z|)sin'"--''(^;. „ sin-''-i(/x, - i|) sin"--''(Ai, + 



3 



where 



^?2(0) 



§1 {x) = ^ (_i)»g"("+i)/2 sin(2„ + d^ix) = q^^^ ^ q"(«+i)/2 cos(2n + \)x. 

Then these equations are written as follows 



dA. nr=iSin''^"'(A. +*|)Bin"-^^(A, -*|) 



na>b sin(Aa - Ab - 

t/.(A„ A.) = -"^-- n"°r^mV^ + '^"^^^ 

3 j — 1 

where fab and yj are given by (O- In Appendix IXI we show that 

UmiXi, ...,Xm) = - 7^)C^i(Ai)C/i(A2)...C/i(A™). 



If we put A, — !■ a;,(/) this is written as 



i-Q/2 

11/ dXjWm(a;i, ... 

7 . J~Ol2 



n 



sm' 



.1, . , 
(xj +i-)sin 



-V] 



"^3 - «o) 



-Ti 0"("-l)/2 n,>, Sin0(xa - - i/afc) 

We investigate the properties of Wi , 

_ ^U0)^2Wx + ^/2)) _ . 2^ . 

^" %92m(0(x + */2)) ~ ttQ'"''^ Q 

where dn(2:, X) is Jacobian elliptic function with 



K{1) 



7r/2 



sin^ fc 



(12) 

(13) 
(14) 

(15) 



(16) 
(17) 

(18) 



This function has poles at x = kQ + i{l — i) with residue (— l)'i and zeros at A = (fc — \)Q + i{l + ^). On the other 
hand X^^-oo ^sech(7r(x + nQ)) has the same poles and zeroes. Then we have 



Using p5|l we have 



Then we have 



Wi(x) = Wi{x + nQ). 



W,„(xi, ...,a;„) = i:)„jWi(xi)Wi(a;2)---Wi(a;,„). 



Wraixi, ...,Xm) ^ ^ Wra{xi + UiQ , X2 + n2Q , ■ ■ ■ 



(19) 



(20) 



(21) 



The weight function in massless case is trigonometric. But in massive case it becomes eUiptic and satisfies the following 
properties, 

"^^jn ('^l 1 ^^2 ; ■ ■ ■ ; -^j ; ■ ■ ■ ; -^fc i ■ ■ M *^m) ~ m (-^1 j -^2 7 ■ ■ ■ 7 -^^A; ; ■ ■ ■ ; i ■ ■ ■ 1 -^m ) J (^^) 

"^'^m ('^l 1 ■ ■ • 1 -^j ! ■ ■ ■ ; •^m ) — '^'^m (^^1 1 ■ • • 1 -^j Q 1 ■ ■ • ; -^m) — ""^^m (-^1 ; ■ ■ ■ ; ~l~ ^1 ■ ■ ■ ; •^m) : (^^) 
"^^m (^1 ; ■■ ■ J -^j 7 ■ ■■ 5 -^m ) — ( 1) ^ ^ ^ ""^^m ( -^l ; ■■ ■ j '^j : ■ ■ ■ j -^m) • (^^) 
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We should note that weight function is doubly periodic but the integrand is not doubly periodic. The correlation 
function has the following properties 



III. REDUCTION TO ONE DIMENSIONAL INTEGRAL AT m = 2 AND 3 

At first we consider EFP defined by This is equal to F(\' ^' '"' ^). As yj = m + 1 — j, fat = 1 and we have 
the integral representation 

m 

l[h"^-'{x,+i/2)h'-'{x,-i/2) 
P{m)^ A'^xW.y,— — , with h{x)=s\n{(j)x)/(j). (27) 



Q/2 YYh{Xa- Xb-l) 



a>b 

Recently there have been an increasing number of researches concerning the properties of the EFP. Particularly, in 
the isotropic limit A ^ l{(f) 0), the general method to evaluate the multiple integral was recently developed by 
Boos et alA .inllll| . 

It is important to note that P{2) and P{3) are related to the nearest and next nearest-neighbor correlation functions 

(5|5^i> = P{2) - 1/4, (28) 
(5|5|+2> = 2(P(3) - P{2) + 1/8). (29) 

We also consider transverse correlation functions (SjSj^i) and (5|5|_|_2)- Then we should evaluate following four 
integrals 

P(2) . r d-xW2^^-\+^/^)^^--//^), (30) 

J -Q/2 h{X2-Xi-l) 



/-Q/2 h{^2 - xi~ i)h{x3 - xi - i)h{x3 - X2 - i) ' 

/•'3/2 h^U + i)h(x2 + ^)h(x2- i^W(x3- I:) , , 

■' J -Q/2 h{X2 - Xi)h(X3 - Xi)h{X3 - X2 - l) 

In Appendix El we transform these to the canonical forms. The two-dimensional integrals of (|3U|I and (|31|l can be 
reduced to one-dimensional ones, 

1 dx 1 



P(2)--+/ (__^+cot((/.x)coth0), (34) 

/oxox \ /-c^+i/a da: 1, X , , cot(0x) , . , 

{S-S^+^)= — . cosh^- . \ 35 

■' ■' J_oo_|_i/2 smhTTX 4 sm'' (px smh0 

These coincide with known results. The main results are that three-dimensional integrals for P(3) and {S^ Sj^2) ca'^ 
also be reduced to one- dimensional integrals. In other words, we have succeeded in performing the integrals for P(3) 
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TABLE I: Estimated values of correlation function (SoSf) from finite rings and our calculation 

iV A = 1.0 A = 1.5 A = 2.0 A = 2.5 A = 3.0 A = 3.5 A = 4.0 

20 0.06135072260 0.1031381904 0.1473943448 0.1795178252 0.1993399825 0.2119381546 0.2204021144 

22 0.06123184969 0.1033007174 0.1480705446 0.1799480453 0.1995078997 0.2119982680 0.2204241750 

24 0.06114201704 0.1034578190 0.1486278313 0.1802452204 0.1996070522 0.2120292690 0.2204342780 

26 0.06107247281 0.1036073511 0.1490859151 0.1804501423 0.1996656632 0.2120452914 0.2204389106 

28 0.06101753236 0.1037485416 0.1494616068 0.1805913409 0.1997003229 0.2120535929 0.2204410300 

30 0.06097337258 0.1038813030 0.1497692404 0.1806886032 0.1997209063 0.2120579113 0.2204420230 

extra. 0.06067894953 0.1061687241 0.1508088848 0.1809161322 0.1999716192 0.2120736609 0.2204679674 

exact 0.06067976996 0.1061446753 0.1511283748 0.1809054601 0.1997511670 0.2120626058 0.2204428632 



TABLE IL Estimated values of correlation function (SoSf ) from finite rings and our calculation 

A =1.5 A = 2.0 A = 2.5 A = 3.0 A = 3.5 A = 4.0 

20 0.04336844424 0.02832886249 0.01874771832 0.01321049421 0.00981186886 0.00757459251 
22 0.04318176477 0.02806630512 0.01860734227 0.01315890868 0.00979394498 0.00756813199 
24 0.04303078324 0.02785949276 0.01851259253 0.01312895682 0.00978483286 0.00756521107 
26 0.04290579437 0.02769531699 0.01844841522 0.01311148673 0.00978017739 0.00756388405 
28 0.04280033153 0.02756430324 0.01840481529 0.01310126617 0.00977778854 0.00756328233 
30 0.04270993496 0.02745937089 0.01837513880 0.01309525214 0.00977655328 0.00756300336 
extra. 0.04170186316 0.02715342577 0.01832502644 0.01302827108 0.00976919713 0.00755456686 
exact 0.04171854326 0.02702161935 0.01831051051 0.01308654379 0.00977523196 0.00756276973 



and (S'JS'J+2): 



P(3) = ^ 



oo+i/2 

-oo+i/2 
oo+i/2 

-oo+i/2 



Ax 1 



sinh TTx 8 



' 3 sinh 



sm 



bx) Vsin {(px) 



5 — cosh 



20^ - cot{(j)x){^ 



3 tanh ( 



i cosh 20 



da; 



sinh Tlx 8 



■(- 

Vsi 



3 sinh 



x) Vsin {(jyx) 



1 — 3 cosh 2(f> + cot{(j}x 



sin^{(t)x) sinh2(/i 
3 cosh(2(/)) tanh((/)) 



(36) 



siv?{(l)x) 



sinh (20). 



Using equations H28II29() we have 



(5*1 5*1+2) - T + 



oo+i/2 

-oo+i/2 
.oo+i/2 

-oo+i/2 



Ax 1 



sinh Tlx 2 

Ax 1 
sinh ■nx 4 



(cot(0a:) coth(0) 



/3sinh~ <p , „ A , I ^^ 

I 2 1 ^ cosh 20 — cot(0a::) I 

Vsin ('(A.t'I / V 



sm 



bx) Vsin {(px) 



3 tanh 
sin^(0a;) 



- 4 coth 20 



(37) 



(38) 



.(39) 



IV. NUMERICAL RESULTS AND 
CONCLUSION 

We calculate numerically the second neighbor corre- 
lations (S'qS'I) and (S'pS'f) for finite rings with N = 
20, 22, 24, 26, 28, 30. In Table ^ and ^ results of ^ 
and (|37|l are compared with numerical results. They co- 
incide very well. 

In the previous paper we treated the second neighbor 
correlations for massless region 1 > A > — 1. In this pa- 
per we have shown that almost the same method can be 
applied for massive region A > 1. This is the generaliza- 
tion of A = 1 result 

(5|5|+2) = ^-^ln2+^C(3). (40) 
The multiple integrals for the correlation functions of the 



XX Z chain at the massive region can be reduced to the 
one-dimensional ones in the case of the next nearest- 
neighbor correlation functions. This property will be 
generalized to the third neighbor or longer distance cor- 
relations. In our future work, we are particularly inter- 
ested in calculating the third-neighbor correlation func- 
tions (5JS'J+3) and (S"; S'J+a) for general A. For A = 1, 
they were recently calculated in [iJl 



(^|5^3) = 3^ - 3 In 2 + f C(3) - f C(3)ln2 

-|C(3)^-^C(5) + fC(5)ln2. (41) 



from the multiple integrals. This should be generalized 
to arbitrary A. 
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FIG. 1: Nearest-neighbor and next nearest neighbor correla- 
tion functions for the XXZ chain. We calculated (SjSj^i), 
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APPENDIX A: ANTISYMMETRIC 
DIFFERENTIAL OPERATOR 

Elliptic theta function satisfies the following quasi dou- 
ble periodicity 



-2iA 



^?i(A). 



(Al) 



We consider 

/(Ai,A2) = (^-^)C/i(Ai)C/i(A2). 

Using ljAl|l and i?2(A) — '!9i(A + tt/2) we can show that 
Ui{X) and /(Ai, A2) are elliptic functions which satisfy 

C/i(A2) = -C/i(A2 + ^(/)) = C/i(A2 + t:), 

/(Al, A2) = -/(Al, A2 + ^(/.) = /(Al, A2 + vr). (A2) 

As /(Al, A2) = — /(A2, Al), / has zeros at Xi + kir + il(f> 
as function of A2. It has poles at kir + i{l — 1/2)0 with 
order 2. Then function / should be factorized as follows 



/(A2) = C(Ai) 



^i(A2-Ai)^?i(A2-a) 
^i(A2+#/2)2 ■ 



(A3) 



Here a is another zero. Using ljAl|l we have /(A2 + i4>) — 
exp(2i(a + Al + i(/)))/(A2). By the condition l|A2(l a must 
be 7r/2 — Ai — i^i. The residue of U\{X) at A — —i4)/2 is 
i. Then we have 



lim ^ (A2 + i^)V(Ai,A2) = -iC/i(Ai). 



By this condition we find 



{d' (0))^ 

^^^'^ = ^^^(^^)^l(Al-K^<^m(Al+^0/2 + V2)' ^^'^ 



and therefore /(Ai, A2) = U2{Xi, A2). Thus we have proven (|20|l for m — 2. 
To prove for to > 3 we assume that it stands for Um-i- Define 

^ m-l g Q 

/„(Ai, A„,) = — TT (— Trr-)Ui{\m)Um-i{\i, •-, A™_i). (A5) 

(to - 1)! ^J-^ dXm dXi 

This function has the following properties: 

/m(Ai, Xm) = -/(Al, .■.,Xm + icfy = /(Al, ...,Xra + Tt). (A6) 



As fm is antisymmetric for exchange of variables we know m — 1 zeroes and m fold pole: 

7n— 1 
1=1 



fm{\n) — C{Xi, A,„_i) — — ■j./r,\m ■ i-^'^) 



a must be tt/2 — J^h^i^ Xi — im(j>/2 by the condition HA6|I . Using 



lim (A„,+z^r/„,(Ai,...,A„) = (-ir-iiC/„,_i(Ai,...,A„,_i), 

AjTi— * — Z0/2 Z 



we can determine C, 



C(Ai, A,„„i) — iJ7„i_i(Ai, A„i_i) ^— J ^ ■ (A8) 

(A, + 1,^/2) + n/2) Yl MXi + ^0/2) 
1=1 1=1 

Thus we show 

Um{x^, ...,Xm) = -— r7 n(ar " ar^^'^^'^^'^^'^-^'^^™)- ^^^^ 

for general m stands if it stands for to — 1. Putting Xj = ipXj we obtain H2()|l . 

APPENDIX B: SOME INTEGRAL FORMULA 

For the integration of canonical forms following formula is important, 

/•Q/2 /•Q/2 

/ dxie^'^^'^'Wiixi + ei) ^ da:ie2™'^^i-^iWi(a;i) 

"'-Q/2 "'-Q/2 

d,,e-"^(^.-.)isech(..0 = ^^P(-2»»'^^i) . (Bl) 
2 2 cosh riffl 



Here we have used ifT^ . Next we consider the integral of ^'^^^'t'^i j (^^2i4>x2 _ g2^0xl^ 

rQ/2 p2inij>xi 
I 

'-Q/2 



/Q/2 g2m02;i 
dxidx2^5-^^^-^5^Wi(a:i + ei)Wi(a;2 + £2) 



.Q/2 /•Q/2 exp 2n - 1 xi - X2 , , ^ 

J -012 J -0/2 2svti(l){x2- Xi) 



= ^i(>^4'(-(2n-l)ti + e2) 



'^Z' f"^^^ eMm2n~l)x,-X2)) ,,, , , . , , 

/ „ . \ rWi(a:i)Wi(x2)da:idx2 

-Q/2 "'-Q/2 2 sm (/)(a;2 - + ei - £2) 



Using ^ we extend the region of integration: 



^ _,g.^(-(2„-i)e.+..) r r £^^P(!M?!i^I)a^3ech;rxisech7rx2dx,d.2. 
J_oo "/-oo 8 sm (/)(a;2 - xi + £1 - £2) 

The term sech7ra;isech7ra;2 can be written as (tanh7rx2 — tanhTrxi)/ sinh7r(x2 — a^i)- Putting y\ = x\ +X2, j/2 = 2:2 — a^i 
we have 



-le 



r . d,2exp(-^<^n,2)) ^ r,^(n-i)...(tanh|(y,+y2)-tanhf(,,-,2))d,,. 

J_oo 16smh7ry2Sm(/)(y2 +£1 - £2) J-00 2 2 
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For n 1 the last integral is estimated as follows: 



An-i)c^yi (tanh + ya) - tanh ^(j/i - y2))dyi 

gi(n-l)0j/i _ 2 



{n - l)(t)i Vcosh^ f (yi + yz) cosh^ - y^) 



dyi 



oo gi(n-l)0(t;i+j/2) _ gi(n- l)</'(j/i - J/2) 



dyi 



(n — l)(t>i 



cosh fyi 



sin(n — l)0y2 
(n - 1)0 



cosh 



4 sin(n — 1)02/2 
sinh(n — 1)0 



For n = 1 this is 4?/2- Then only integral with respect to 7/2 remains. To circumvent the poles on the real axis we 
change the path of integral from [—00, 00] to [—00 + i/2, 00 + i/2]. Thus we have 



Q/2 g2zn0a:i 

da;idx2 ^-7 — 

'Q/2 ^ ^ 



= < 



'(2n-l)ci+e2) 



Wi(a;i + ei)Wi(x2 + 62) 

da;exp(-m0x) 
-oo+i/2 4 sinhTra sin 0(a; + ei — £2) sinh(n — 1) 



sin(n — 1)02: 



for n 7^ 1, 



_,Vi0(-ei+e2) 



oo + j/2 

-oo+i/2 4sinh7ra;sin0(x + ei — £2) 



(B2) 



da; cxp(— i0a;) 



a;, for n — 1. 



If the integrand is given by a canonical form we can estimate the integral by these formulae. 



APPENDIX C: CANONICAL FORM OF INTEGRAND 



The weight function W,„ is elliptic function with periodicity Q,2i and the integrand is periodic function with 
periodicity Q. Then integrand is expressed by zi = exp{2i(j)xi +0). If integrand has the form (zi — l)™/(zi, Z2, Zm) 
and / has no singularity as a function of zi, the m-th pole of W™ at xi = i/2 is cancelled and shift of path from 
[-(3/2, Q/2] to [-Q/2 + i, Q/2 + i] becomes possible: 

Q/2 rQ/2 

dxiW„(zi - l)"7(^i, ^2, ...) = - / dxiW„(g2zi _ l)™/(g'^i, Z2, ...)• (CI) 

-Q/2 "'-Q/2 

Here we have used (|23|l . In the same way we have 

/•Q/2 /■Q/2 

/ d.TiW„(<z2zi-l)™/(zi,Z2,...) = - / dxiW„(zi-l)"7(g-2^i,Z2,...)- (C2) 

J-Q/2 J-Q/2 

Actually / might have single pole at zi — q^"z2, n = integer, because Wm has single zeros at these points. Then 
we have 



Q/2 C _ l\ra j.Q/2 , 2 _ 

da;iW™ii ^/(^,,^2,...) = - / Ax^^J^ >f(q^z^^z^^,„)^ (C3) 

-Q/2 Zi - q^Z2 J-Q/2 q [Zl - Z2) 

Q/2 (- 2 -|\m ^•Q/2 / -|\rn 

da;iW„^V ^—f{zi,Z2, ...) = - / dxiW,„ii )-f(q-^^,^,,^ ...). (C4) 

-Q/2 q Zi - Z2 J-Q/2 Zi - Z2 



From (j221l we have 



Q/2 /•Q/2 

d™a;W,„.g(zi,Z2,Z3,---) = 

-Q/2 J-Q/2 



/•Q/2 

/ d™a;W™g(z2,zi,Z3,...). (C5) 

"'-Q/2 



Suffix 1 and 2 can be replaced by general j ^ k at 1 < j , k < m. These rules are sufficient to simplify the integrand. 
Hereafter we write A B ii 

Q/2 rQ/2 

d™xW„A= / d"'xW^B. 

-Q/2 J-Q/2 
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At first we consider two-dimensional integral. Equation H30() and H31(l are rewritten in terms of q = exp(— </>) and 
zi = exp{2i(l)xi + 4>), 



P(2) 

The integrand is transformed as follows 

{q^zi - 1)(Z2 - 1) _ 1 



da;idx2W2— — a"^' 



, (-^2-1)% 1 , ((?2z2-l)\ 



2i(f){z2 — q^ zi) 2i(() Z2 — q^zi' 2i(f> 

We can show that z'l/{z2 — zi) ^ — zi + 9^^/(22 — ^i), because 



g2(z2 - zi) ■ 



(zi - 1)2 {q^z, - If 



Z2 - zi Zi- q-'zi 



q^zi + Z2~2- 



(^2 - 1)^ 
Z2 - q^zi 

2, 1 \2 



-{l-q')zi + 



2, , (g"z2-l)^ ,^ 2^ , (g'^l-l) 



^2(22 - Zi) 



-(l-<Z^)zi + 



g2(z2 - Zi) ' 



Then we have the canonical form 



2l0 Z2 — Zl 

Using ||BlJ and we evaluate /^^/2 dxidx2Wi(xi + ei)Wi(x2 + e-,) ^^ : 



-2i4>ti 2 /•00+V2 



da; 



4z(/) 40 J-oo+i/2 sinh7ra;sin0(x + ei — £2) 
Using H2U|I we have 



^^-i^(^+e,-e2) _^ gi0(ei+62) coth^sin^a; 



p(2) = lim (^2 " aoinini = dMii- 



(|31|) is rewritten 



The canonical form of integrand is 



"^^'h H w (g^^l-l)(^2-l) 
dxida;2W2 — -—7 r — . 

-Q/2 Al(j)q[z2-Zi) 



1 -(l + g2)zi+2 



4z(/)g Z2 — Zl 



Then we have 



(5J5J+i)= hm (92-91)- 



oo+i/2 



(C6) 



(C7) 



(C8) 



(C9) 



dx 



-00+1/2 sinh7rxsin(/)(a; + ei — £2) 



xcosh (/>e 



-i0(a:+ei-e2) _ „j0(ei+C2); 



sinh0 



Equation is 



Q/2 



(9'^1-1)'(9'^2-1)(Z2-1)(Z3-1)^ 



/-Q/2 (Z2 - g2zi)(z3 - gf2zi)(z3 - g2^2) 

The integrand is devided as follows 

{q'z, - l)\q'z2 - 1)(Z2 - 1)(Z3 - 1)2 _ (q^z, - ifiq^z^ - l)(z2 - l)(z3 - 1)^ 



(CIO) 
(Cll) 

(C12) 

(CIS) 

(C14) 
(C15) 



(Z2 - q^zi){z'i - g2zi)(z3 - q^zi) 
/I 1 



Z2(l-92) 



1 



■ (z3 - q^zi){z-i - q^zi) (Z3 - q^Zi){z2 - q^Zi) (Z3 - q'^Z2){z2 - q^zi) 
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The first and second terms are symmetrized by the exchange of variables and the third term is devided to four terms 



(Z1-1)(Z2-1) 



{q^zi - 1)^(1 + q^Z2Z3 ~ Z2- z:i) , {q^ Zi - 1)^Z3 - if 



Z2Z'i{l - q^){z2 - q^ zi) Z2{1 - q^) 

{q^z,-lf{z:,-lf , (Z3-I)2(q2^i-1)3 {z:^ ~ if {q^ z, ^ if 



Z2{1 - q'^){z3 - q^ Z2) Z2{l-q'^){z2-q'^zi) Z2{1 - q'^){zs, - q^ Z2){z2 - q^ zf}' 
This is equivalent to 

' [q^z, - 1)2(^3 _ 1)2 (q^z, - 1)2(Z3 - 1)2 _ [q^z^ - lf{z, - ifq^ 



(C16) 



^2(l-g2) q^z^Z2{l-q^) 
jq^z, - l)(zi - 1)((Z3 - 1)^ - jq^z, ~ if) ^ {q^z^ - Ifzs - 1) {q^z, - if - (z^ - 1)^ 
Zi(l-gf2) 23(1-^2) Z2-q2zi 

(g2zi _ 1)3(^3 _ 1)3 (Z3 - l)3(g2^i _ 1)3 (^^ _ i)3 (^2^^ _ i)3 



q^ziZ2{l - g2)(z3 - q'^z2) Z2Zj,{l - q'^){z2 - q'^zi) Z2{1 ~ q'^){z3 ~ q^Z2){z2 - q^zf 
We use following formula to simplify terms in the second parenthesis; 



(C17) 



(^3 - 1)^ . {q^Z^-lf (^2^)3 (zi-1)^ . 
2—f{zl, Z2) 2? ^2 ' 2 /^2, ^3 27 \ /^2, ^3 ■ C18 

Z3-g22;2 92(23-2:2) z2-g2z;^ 92(22-21) 

The first parenthesis is simplified only by the exchange of variables. Then we have 

'3g2z223 z2 (^2^^„^)3(^2^^_1)3_^4(^^_1)3(^^_1)3 



jni3^f^^-3(l+9^)^+8^ 

V z-} zo Z-) 



Z2 Z2 Z2) V g4(i _ ^2)2223(22 - Zi) 

(q2z3 - l)3(zi - 1)3 



(C19) 



Z2(f{l - q'^){z3 - Z2){Z2 - Zi), 

The last term in the second parenthesis is 

(g223 - 1)3(2, - 1)3 _ (g2^^ _ 1)3(^3 _ i)3 
(72(1 - 52)22(2:2 - Zl) 

The terms of the form /(zi, 23)7(22(22 — 21)) can be replaced by — /(zi, Z3)/(ziZ2) + /(zi, 23)7(21(22 — zi)). Then we 
have 

m - f ^ ~ 3(1 + ,2)fi + 3fl) _ i,^z,-mq^z^-lf-qHz3-mz,-lf 



Z2 Z2 Z2I g"*(l - g2)2iZ223 

(q223 - 1)2(21 - 1)2 + (q2^3 _ i)(^3 _ I)(g2^^ _ i)(^^ _ 1) + (^2^^ „ 1)2(^3 „ 1)2 



g22i22 



Z2 - Zl 



+ + + 1)^1^3 - 3(g* + g' + 1)^1^3(21 + 23) + 3(g2 + i)(22 + 32123 + 22) 



^1 3,1 1, l + g-4 

(— + — + 9 21 + 23 + ^ — + — — 

Z3 zi q^ 21 23 2123 

(g223 - 1)2(21 - 1)2 + (g223 - 1)(23 - l)(g22, - l)(2i - 1) + (g22, - 1)2(23 - 1)^ 

g2zi 



(C20) 



Then all terms are in the form of z1z\z^ or 2123/ (22 — 21) with 2 > a, &, c > — 1. We should note the following relations 

2f( ^ ^ 3(1 + g2) _6 ^ l + g2 

Z^]{Z2,Z3) - I ■ 21 - + )/ 22,23 , 

\ 1 + g^ 1 + g^ g'^(l + g*)2i 

^/(^3)~(^^^ r- 4/^^\ -^i)/(^3). (C21) 

22-21 Vg2(22-2i) g4(22-2i)2i / 
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By these relations we can put powers a, 6, c are at 1, 0, —1. The canonical form is 
{l + q^zi 1 rl2(l + g2) {l + lOq^ + q'^)zi 3(1 + 



q^ Zj, {Z2 - Zl) 



, l + lOg^+g^ , , 2 3(1+<Z2)| 1 I 3(l + g2)2 3(l + g2)zi (l + g2)(l+g2^^4) 

+ ^3 < 2 1 + 9 )z\ + 



q^zi J Z3 I 9'' <f <fzi 

(C22) 

Multiple integral can be done through (|B1IIB2|) : 

—3 / d3xWi(xi + 6i)Wi(a:2 + e2)Wi(a;3 + ea)^ = » + -3 / — —— r 

H J_Q/2 IH ^(P J-oo+i/2 Sinh7r2;sin0(x + ei - £2) 

+ .,<^(e,+£,) Sin'/'a: _ (1 + lOg' + g^) ,^(_^_,, 3(1 + g^)^ .z<»fa+3.i +e,) sin 20a 



e 



g sinh g2 g2 sinh 2(/) 

^2^0.3. 1 + lQq^ + qi sin^x , 3(1 + g^) , , . , +^ , 3(1 + g^) , ^ . sin 2</)a; 



r 1 + +g ^»^(,,+e,) Sm0x ^ £ii_+gj^,^(_,_,^+,,)^ ^ dil_+gJ^,^(,+3,^+,^)Sin^| 
I sinh (h Q a sinh 2(i) J 



cosh I g2 sinh (p q q sinh 20 . 

^ r_3(l + g^^,^(,^+,^)Sin0x ^ 3(1 + ^ ^.^(-.-e.+e.)^ ^ (1 + g^)(l + g^ + g^) ^,^(^^3^^ ^^^^ sin 20x | 

cosh g2 sinh (j) q q^ sinh 20 J 

(C23) 

Then we have 

P(3) = i lim (a2-9i)(53-ai)(93-52)in2Sl = dSni. (C24) 

2 ei,e2,e3^0 



Equation (|33|l is 



'-Q/2 g^(^;2 - ^l)(^3 - Zi)(z3 - 9222) 

The integrand is devided as follows 

(g2zi-l)2(g2z2-l)(z2-l)(23-l)V 1 



-(-1 + 



g2(l - g2)2:2(z3 - g2z2)(z2 - Zl) g^(l - g^)z2 (^2 " ^l) ^3 - 9^22 ' 

(g^Zl-l)^(z2-l)(^3-l)^ , (g^Zl - 1)^(Z2 - l)(g^^3 - 1)^ 
g2(l - g2)z2(z2 - Zl) g''(l - g^)z2(z2 - Zi)(z3 - Z2) ' 

Using the equivalence 

/(2;i,2:2,Z3) 1 //(zi,Z2,Z3) - /(zi,Z3,Z2) 



(Zl - Z2)(zi - Z3) Zl - Z2 V Z2 - Z3 

we replace the term like (zi — Z2)(zi — Z3) in the denominator by zi — Z2. The canonical form is 



261 



Z2 - Zl 



(l + g2)(3 + g4) (3 - 2g2 + 3g4)zi (l + g2)2 



g* g2 g^Zi 



(C26) 



g2(l-g2)2;2 V (23 - Zi)(z3 - g2z2) (z3 - Zl ) (z2 - Zl ) (z3 - g2z2) (2:2 - Zl) 

The first term is transformed as 

(g^Zi-l)^(g2z2-l)(z2-l)(z3-l)^ ^ (g^Zi-l)^(z2-l)(g^Z2-l) ^ ^2^^ _ 3 + (^'^^ " ) 
g2(l - g2)z2(z3 - Zi)(z3 - g2z2) g^ (1 - g2)z2 (z3 - ^1 ) Z3 - g2z2 

(g^Zi-l)^(z2-l)(g^Z2-l) , 2 (g^Zi-l)^(g-^Z2-l)(z2-l)^ 

^ 271 2^ — 7 ^ l^s + 9 ^2 - jz — 7 rj 

g2(l - g^)Z2(Z3 - Zl) (1 - g^)Z2(Z3 - Zl)(Z3 - Z2) 

The third term is transformed as 

(g2zi-l)2(g2z2-l)(z2-l)(z3-l)2 _ (g2zi-l)2(z2-l)(z3-l)2 Z3 - 1 , 



(C28) 



(C29) 



;6-2g2 (l + g^)zi , (l + g^)(-2 + g^) ] , 1 f 4 + 4g2 2, l + 9'\ 

^3 i —3 75 + TaZ r + T 1 (1 + 9 )^i - — I— } ■ (C30) 
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Then the integral of this form is 

i [Ql'^ 

— ^ / d^xWiCa;! + ei)Wi(a;2 + e2)Wi(a;3 + e^)^^ 

^4>^ J-Q/2 

<.oo+i/2 



1 



dx 



-oo+i/2 

2\fi I ^4 



sinh7rxsin0(x + ei — £2) 



cosh( 

cosh(/ 
Then we have 



(l + g-')(3 + g'^) ^,,^(,^^,^) Sin0x ^ (3_2g^_F3Og,0(-x-ei+.2)2. + (1 + 9 )\ .»0fa+3ei+e.) sin 2,?!)a: 
(7-^ sinh (/) q"^ sinh 20 

g'''^^^ /6-V^.^(e,+..)£ni0x _ + , (l + g^)(-2 + g^) ^,,(,^3,^^,,^ sin20a; 



sinh (j) L 

g sinh0 5 g sinh2(/) 



sinh 20 



(C31) 



(C32) 
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